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Abstract 

The problem of pairing in anisotropic electron systems possessing patches of fermion 
condensate in the vicinity of the van Hove points is analyzed. Attention is directed 
to opportunities for the occurrence of non-BCS pairing correlations between the states 
belonging to the fermion condensate. It is shown that the physical emergence of such 
pairing correlations would drastically alter the behavior of the single-particle Green 
function, the canonical pole of Fermi-liquid theory being replaced by a branch point. 



PACS: 71.10.Hf, 71.27.+a, 74.20.Mn 

The ground state of conventional superconductors at T = is known to be a condensate of 
Cooper pairs with total momentum P = 0. In Fermi-liquid theory, the familiar BCS structure 
of the ground state is associated with the logarithmic divergence of the particle-particle 
propagator at P = and is independent of the details of the pairing interaction. However, a 
markedly different situation can ensue in strongly correlated systems in which the necessary 
stability condition for the Landau state is violated and the Landau quasiparticle momentum 
distribution suffers a rearrangement. Under certain conditions this rearrangement leads to 
a fermion condensate (FC) - a continuum of dispersionless single-particle (sp) states whose 
energy e(p) coincides with the chemical potential fi over a finite (and in general disconnected) 
domain p e Q in momentum space [1-11]. As a result, the preference for pairing with P = 
comes into question because of the degeneracy of the FC sp spectrum. In this case, the 
nature of pairing depends on the configuration assumed by the FC. 

Here we study a two-dimensional square-lattice system in which the FC is situated in 
domains adjacent to the van Hove points, while the sp states with ordinary dispersion are 
concentrated around diagonals of the Brillouin zone || pj[ . To begin with, we focus on 
the nature of particle-particle correlations in the FC subsystem and ignore contributions 
from the sp states with nonzero dispersion. Traditional BCS singlet pairing correlates only 
the sp states belonging to diagonally opposite patches of the FC; the description therefore 
involves the single collective operator C p = a p ._a- p , + and its adjoint CL which connect the 
ground state with states of N and N =p 2 particles. However, in anisotropic electron systems 
inhabiting crystalline materials manifesting fermion condensation, all four FC patches should 
be treated on an equal footing. Hence an additional relevant collective operator Q p = 
a P)+ a -p+Q)- ; enters the picture, together with its adjoint Q p . With Q = (tt/1,tt/1) where I 



1 



is the lattice constant, this operator characterizes the pairing correlations affecting sp states 
located in the neighboring FC patches. If such additional correlations are involved together 
with the ordinary BCS correlations, the ground-state wave function evidently loses the simple 
BCS structure. 

A salient feature of this nonabelian exemplar of the pairing problem is the presence of 
two degenerate collective modes in the particle-particle channel. The creation of a C-pair 
(Cooper pair) by the operator C p , followed by subsequent annihilation of a Q-pair by the 
operator Q p , gives rise to an excited two-particle-two- hole state of the iV-fermion system. 
As we shall demonstrate, this process enmeshes a whole band of many-particle-many-hole 
states and changes the structure of the single-particle Green function dramatically. 

We restrict considerations to the simplest, <5-like form of the interaction in the particle- 
particle channel, with strength parameter A. Also, we assume that all the particle-hole 
contributions have already been taken into account in terms of an effective single-particle 
Hamiltonian having sp spectrum e(p) Accordingly, only pairing contributions should be 
incorporated in the equation for the Green function G a/3 (x,x') = -i (T^(x)^(x')). This 
equation, derived with the aid of equation of motion [e — e(p)] ip a (x) — X^(x)^j(x)tp a (x) = 
0, takes the form 

(e - e(p))G aft (x,x') + iX(0\T^(x)^(x^ a (x)^ p (x)\0) = 5 ( x ~ x ') ■ W 

In the ordinary pairing problem, the average (0\T^(x)^(x)i(} a (x)^(x')\0) is decoupled as 
(0\Ta pa (t)a-^ a (t)\C}(C\Tal uy (t)al pul3 (t + T)\0}. In the generalized case being developed, 
where the ground state is connected with two different states given the labels C and Q, this 
same average has the extended decomposition 

(0|Ta Pl , 7 (t)a_ Pl , a (t)|C)(C|TaJ )i7 (t)a t _ P)/3 (t + r)|0) + 

+ (0|Ta Pl , 7 (t)a_ Pl+Q , a (t)|Q)(g|Ta pi7 (t)a t _ p+Qi/3 (t + r)|0). 

For simplicity we henceforth omit spin indices a, (3, 7, etc. The equation for the Green 
function then reads 

G(p,e) = G o (p,e)[(l-AF+(p,e)- J DF +(p, e )] , (2) 

where G Q (p,e) = [e — e(p)] -1 . In the time domain, the quantity F^ has the expres- 
sion F^ Q (p,t) = (C\Ta p (t)aL p (t + t)\0) and is interpreted as the transition amplitude 
between the ground state and a state differing from it by the presence of a single C- 
pair. Similarly, F^^p, r) = (Q\Ta^(t)a + n(t + t)\0) is the transition amplitude be- 
tween the ground state and a state differing from it by a single Q-pair. The diagram- 
block A oc (0\Ta Pul (t)a- Pl!a (t)\C) oc Fi t0 (r = 0) has the same meaning as the gap or- 
der parameter of BCS theory, while the new ingredient D has the corresponding structure 
D oc (0|Ta pi)7 (t)a_p 1+ Q ia (f)|Q) oc F 0>1 (r = 0). 

Employing the complementary equation of motion [e + e(p)] ^\ t (x) + \'ip\ l (x)%l)\ 1 {x)'ip 1 (x) = 
0, one can derive equations for the transition amplitudes F{|" and F^: 

[e + e(p)] (C\Ta p (t)al p (t + r)\0) = -A£ [(C|ra pl (t)al pi (t)|0)(0|Ta p (t)a p (t + r)\0) + 

pi 
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+ (C|Ta+(t)al pi+Q (t)|Cg- 1 )(Cg- 1 |Ta p (t)a+ +Q (t + r)\0)\ , (3) 
and analogously 

[e + e(p)] (Q\Tal(t)al p+Q (t+r)\0) = -A£ [(Q\Tal(t)al pi+Q (t)\0)(0\Ta p (t)al (t+r)\0) 



pi 



+ (Q|Tat i (t)a t _ p+1 (t)|gC- 1 )(gC- 1 |Ta p (t)a t _ p+Q (t + r)|0) 
Equations (^) and (^) are conveniently rewritten in the compact form 



(4) 



F + 
F + 



G, 



A + G + D + Fi-i 
D + G + A+F_! ! 



(5) 



where G Q = — [e + e(p)] . In contrast to BCS theory, the above system of equations 
for the key quantities G, -Fi.o, and Fq.i is not closed: new transition amplitudes iq,_i ~ 
(CQ-^Tap^alp^t + r)\0) and F_ M ~ {QC- l \Ta Pl {t)al ?1+Q {t + r)\0) come into 
play, for which further equations must be derived, and so on ad infinitum. 

Upon neglecting the contributions from the new amplitudes and F-i t ±, the usual 

BCS equations are recovered. If we keep these ingredients, the situation is complicated but 
still amenable to analysis. In the next step, we form equations for the amplitudes and 
i by repeating the process that led to Eqs. (0) and (|5]). Ignoring the variation of the 
blocks A and D due changes of the states between which the transitions occur, we have 



iVa = -G o [DF+ + AF+_ 1 

F-1,1 = -g \af++df+ 



(6) 



Suppose now we neglect the last term on the right in each of these equations. Then the 
system of equations (Q), @, and (Q) closes, and may be solved without difficulty In this 
approximation, one finds 



G -AG( 2 7(1 + K Q ) , (F$)+ = G-DG^/(1 + K c ) , 
where 

K c (p, e) = G (p, e)A+G-(p, e) A , K Q (p, e) = Gj(p, e)DG (p + Q, e)D^ 
Assembling the Green function G via Eq. (Q), we obtain the approximant 

G n 



(7) 



G^ 



with denominator 



= 1 + 



K, 



c 



+ 



1 + Kn 1 + K, 



c 



(9) 
(10) 
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Equations for the transition amplitudes F2-1 and F_i 2 , omitted in the foregoing manip- 
ulations, may be derived along the same lines with the results 



F, 

F + 
r -1.2 



2,-1 
+ 



G (j 

G: 



D+F„ lA + A + F_ 2 , 2 



Again deleting the last term on the right of each equation and proceeding as before, one 
finds 

G-DcG&> <3v + G-D Q GW 

1^0,1 J _ 



(3)> 



1 + 



and 



with 



l+K c 



! + I 



_G_0_ 

Z( 3 ) 



= 1 + 



K, 



c 



1 



+ 



Kr 



1 + 



(12) 

(13) 
(14) 



1+K C * ' 1+-KQ 

Indefinite continuation of this process generates the coupled set of equations 



F 

± n,—n 

F 

± —n,n 

7+ 
n,—n+l 

7+ 
— n,n+l 



-G AF n+1 + DF n n+1 



-G 



AF 



— n+l,n 



n,n+l 



G 

G: 



D + F_ r 



A + -F n _ 1) _ n+1 
A + -F_ n _ ljn+1 



etc. 



(15) 



Together with Eq. (0) for the Green function G, this system can be written in closed form 
in terms of two new Green functions Gc and Gq\ 



G 

G c 
Gn 



G 
G 
G„ 



G A+G C AG - G D + G Q DG . 
G D + GqDGc , 
G A+G C AG Q . 



The system ( p!|) has the solution 

G 



G 



where 



S + K Q /S 



5 = 1 + 



G"AG 



G A+G -A 



G-PG D+ 

1+G A+G- A... 



Kq/S 



1 + 



Foa 



G-DG 
S 



(16) 



(17) 



K, 



c 



1 + K Q /S ■ 



The latter equation, rewritten as S 2 + (Kg — K c — 1)5 — Kq = 0, is solved by 



S 



Kc-Kn + 1 



-^K c -K Q + lf + AK Q 



(19) 
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the plus sign of the square root being chosen so to give S = 1 if Kc = Kq = 0. With the 
evaluation 



Kn 2K 



Q 



S K C -K Q + 1 + ^(K c -K Q + \y + AK Q 

1 

~2 



K c - K Q + 1 - J(K C -K Q + l) 2 + AK Q (20) 



we have 

1 + Kq/S = (K Q -K c + I) /2 + yJ(K Q -K c + l) 2 /4 + K c . (21) 

It may be noted that the factors S and 1 + Kq/S entering the formulas for F 01 and F 10 
transform into each other under interchange of Kq and Kq. 

Inserting relations ([19]) and (f2(|) into Eq. flT7|), we arrive finally at 

G = 7777 77 . .xn . . t f 1/2 = 7 7775" (^) 



[(K C -K Q + 1) 2 + 4K ( 



if 2 + K> - 2/^q + 2K C + 2^ Q 



1/2 



and 



'1,0 



G n AG G n DG 



(1 + Kq/S) [(K c -K q + l) 2 + 4K Q } 1/2 ' ' S \(K C - K Q + l) 2 + AK Q } 1 ' 2 ' 

(23) 

The conventional Fermi-liquid-theory pole in the Green function is no longer present; it 
has been replaced by a branch point. This is our primary result. The Green function G(e) 
possesses a nonzero imaginary part over a finite interval in e delimited by the two zeros of the 

Too "I 1/2 

denominator Z = K c + Kq — 2KqKq + 2Kq + 2Kq + 1 .In fact, this result is generic 
in the sense that it is independent of the details of the interaction in the particle-particle 
channel and certainly not specific to the assumed 5-function form. Such details only affect 
the form of the equations for the gap functions. 

Beyond this result, there is the obvious question of whether non-BCS pairing of the kind 
described here can win the contest with ordinary BCS pairing. To decide this issue, one needs 
to know the gap functions A and D. For our simple model interaction characterized by a 
single parameter A < representing the coupling strength, these functions obey equations 
of the customary form 

A = -A / Fi, (p,e)^p— , D — —X F 0A (p,e)dv p — . (24) 
J 2m J 2m 

Here dv p denotes the usual momentum-space volume element, and at nonzero temperature 
T the factor tanh(e/2T) is to be inserted in the integrands. In more realistic models, the 
constant A is replaced by appropriate blocks Vi 5 o and Vo,i of Feynman diagrams irreducible in 
the particle-particle channel, and the behavior of these interactions in position or momentum 
space (e.g. with strongly repulsive as well as attractive components) may lead to different 
conclusions than those drawn below. 
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The solution of the gap equations Q2H]) requires a knowledge of the sp spectrum £(p) = 
e(p) — \i. As a result of pairing correlations, e(p) differs from zero even in the momentum 
region occupied by the fermion condensate @ . The truth of this assertion becomes transpar- 
ent when the BCS case specified by D = and A = Ac is considered. The first of Eqs. (F24|) 
then reduces to 

Ac = -A J ^7^dv p = -\J Jn(p)(l-n(p))dv p . (25) 

We proceed under the assumption that pairing does not significantly affect the momen- 
tum distribution n Q (p) which has been determined for the "normal" system through the 
variational condition 5Eo/5n(p) = fi, the presence of a FC having been established. Re- 
placing n(p) by no(p) in Eq. (|25D , one obtains the relation n (p) = [E(p) — £(p)]/2£7(p), 
where E 2 (p) = £ 2 (p) + A 2 C . It follows that E(p) = A c / [4n (p)(l - n (p))} 1/2 and £(p) = 
(1 - 2n (p))A/ [4n (p)(l - n (p))] 1/2 §. Thus, inclusion of pairing correlations results in 
an inclination of the FC plateau in the spectrum £(p). 

Next we argue that within the current model involving a single interaction parameter 
A, the energy gain produced by BCS pairing exceeds that coming from the new pairing 
scenario. To facilitate the demonstration, we neglect the insignificant difference between 
£(p) = e(p) — fi and £(p + Q) = e(p + Q) — /i in the FC region. One can then verify that 
the system of equations (|17D and (|24]) has the solution A = D = A N . For this case the first 
formula in Eq. (|17]) simplifies to 

G = -. Go , (26) 

and we are led to 
G(p,e) 



£ + £(p) tc\ f £ + £ de 



\[e- £(p)][e 2 - ^(P) - 4A^] ' ^ J J e V [e - + ^% ~ e 2 } * ' 

(27) 



and 

-5 



A, = A / 2An ^dv p , (28) 



-E 



where E(p) = [£ 2 (p) +4A 2 ]. 

To determine which of the two competitors - BCS or non-BCS - prevails, we should 
compare the superfluid corrections SE^ and 8Eq to the corresponding ground-state energies. 
For this we employ 

r\ A 2 4- D 2 M A 2 

5E N = - ?±-±£-d\'=-2 £%-d\', (29) 



Jo (A') 2 Jo (A') 2 ' 

a result derived in the same manner as the analogous formula for 8Eq appearing in the 
Landau-Lifshitz textbook |12|. Upon further analysis it is found that at T = the gap 



parameter A at is markedly smaller than Ac, which decisively favors the BCS solution. 
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However, the energetic advantage of the conventional BCS solution decreases as the 
temperature increases, and it disappears entirely as the critical temperature T c is approached. 
Consequently, if the coupling constant for the states in the neighboring FC patches were to 
exceed that for the states in the opposite patches, one could encounter the same situation 
as in liquid 3 H: close to T c , superfluid 3 He is in the A-phase, but upon cooling the liquid, 
the A-phase surrenders to the B-phase. Even so, if contributions from sp states with normal 
dispersion are taken into account, the conditions for the ascendency of non-BCS pairing 
become more severe, suggesting that such exotic phases are very rare visitors. 

We are indebted to M. V. Zverev for numerous discussions. This research was supported 
in part by NSF Grant PHY-9900713, by the McDonnell Center for the Space Sciences, and 
by the Russian Fund for Fundamental Research, Grant No. 00-15-96590 (VAK). 
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